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INTRODUCTION 


The interaction of gases with solid surfaces is a problem of considerable 
practical importance and has received much attention in the literature. Application 
of the solution of some of these problems is found in surface coatings, vacuum tech- 
nology, thin-film growth, and surface contamination. In addition, new and interest- 
ing physical processes are being studied through surface scattering experiments. 

Our principal interest is in the study of the growth of thin, solid films on salt 
surfaces . 

The interaction of gases with solid surfaces is characterized by a long-range 
attractive force combined with a short-range repulsive force (ref. 1) . It is usually 
assumed that the total gas-surface interaction can be approximated as the sum of the 
gas atoms and the individual atomic constituents of the solid over the two-body 
interaction potentials (ref. 2). Thus, determining the interatomic potential is of 
fundamental importance . 

In examining the literature on the binding potentials of individual atoms on 
salt surfaces , one is amazed by the diversity of approaches and by the wide range of 
final values obtained by various authors. We found no. consistency in the treatment 
of the long-range attractive forces even though these forces are reasonably well 
defined from atomic spectral properties . The short-range repulsive core parameters 
were taken from diverse studies on gas properties and on ionic crystals with no 
attempt to treat all the data on a systematic basis. After a brief review of the 
theoretical developments in this field of research, we present a program to con- 
sistently analyze and develop a data base of physical interaction parameters for 
future studies of the interaction of noble metals with salt surfaces. 


REVIEW OF THEORETICAL DEVELOPMENTS 

Models of the interatomic potential are usually approximated with a long-range 
attractive-force potential varying as the inverse of the separation of the atoms 
raised to the sixth power. The attractive potential is usually identified with the 
London dispersion force (ref. 3) . The short-range repulsive core potential in early 
models was taken as the inverse of the separation raised to the twelfth power by 
Lennard-Jones (refs. 1 and 2) or as an exponential by Buckingham (ref. 4) . The 
Lennard-Jones potential 


<l>Lj( r ) 



has two free parameters usually taken as the equilibrium well depth £ and the mean 
equilibrium separation R, where r is the separation of the atoms. (A list of 
symbols appears after the references.) The zero potential crossing radius 0 is 
fixed at 0.89R, independent of well depth. The constant of the London dispersion 
force y is then given by 2eR^. The Buckingham potential is usually written as 




where 3 is a parameter characterizing the repulsive force. The corresponding 
value of the constant of the London dispersion force is given by £R 6 /[1 - (6/3)], 
which is equal to the Lennard- Jones value only for 3 = 12 . The zero potential 
crossing radius is 0.89R when 3 ~ 14.5. Since the second virial coefficient is 
related to O 3 , R 3 , and £, most gases yield O x 0.89R and result in a London dis- 

persion constant, which is 18 percent lower for the Buckingham potential than that 
found for the Lennard- Jones potential. The analysis of gas thermodynamic data shows 
more consistent results from use of the Buckingham potential (ref. 5) . Early work 
on interatomic potentials was largely devoted to extracting potential parameters 
from experimental data and searching for semiempirical relations among parameters to 
provide some semiempirical predictive capability. 

The long-range attractive forces between atoms were first analyzed by London 
(ref. 3) . He showed that the first-order perturbation between two atoms was zero 
but that second-order perturbation theory gave rise to an interatomic potential 
related to the atomic multipole transition amplitudes, giving rise to the name dis- 
persion forces. The r -6 term was related to the interaction of an induced dipole 
in each of the two atoms and the potential coefficient was related to the dynamic 
polarizabilities. He further showed that if the dipole polarizability was replaced 
by a single-pole term at frequency w 1 , then the dispersion force constants could be 
written as 
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where and are the corresponding atomic static polarizabilities of the two 

atoms. (The unprimed and primed symbols indicate the two atoms.) Since the energy 
poles appear near the ionization potentials I, we have the approximation 
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Another approximation used in many surface physics calculations is that of Kirkwood 
and Muller (refs. 6 and 7): 


KM 


6mc 


a 2 a 2 


(a 2 /x m ) 


( a^/x;) 


where m is the electron mass, c is the velocity of light, and X m and X^ are 
the magnetic susceptibilities of the two atoms. These formulas may be used to derive 
general relations between the force constants ]i of interaction potentials for dif- 
fering atomic systems (ref. 8), denoted by subscripts 1 and 2, as follows: 
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which can be approximated as 


P 12 ~ \/ y ll y 22 


It can be seen from this last equation that 
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which follows from the Lennard- Jones potential and from the fact that R-^ ~ \/ r ii r 22’ 
The combination rules for the Lennard- Jones potential are the same as for e q 2 ' and 
consequently, we can derive 
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The exponential repulsive term was proposed by Zener (ref. 9) , who relates the 
exponential decay constant to the asymptotic electronic wave functions 


ip = exp^-^2I/a Q e 2 x 

where I is the ionization potential, a Q is Bohr's radius, x is the electron 
radial coordinate, and e is the electron charge. Zener's results suggest that 


g 12 = ifhl + ^22 

R 12 2 \ R 11 R 22 


is an appropriate combination for the exponential parameter 3 of the Buckingham 
potential . 

The above formalisms have been used to study physical adsorption on surfaces 
over the years. Studies of noble gas adsorption on salt in particular (ref. 10) 
relied on salt cohesive- force data (ref. 11) . The London dispersion force was 
usually approximated through use of the Kirkwood-Muller formula, the magnetic sus- 
ceptibility data of Hoare and Brindley (ref. 12) , and the polarizability of Mayer 
(ref. 13) . The repulsive core parameters were estimated from virial data for argon 
and studies of salt crystals (ref. 11). In addition to these forces, the electro- 
static attraction caused by the induced dipole moments in the adatom was added 
(ref. 10) . Note that static polarizability was estimated from the interatomic 
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potential data and may or may not reflect the dispersive properties of isolated 
atoms. Later work used the same methodology with improved data (ref. 14). 

In most applications of surface physics, the above historical development is 
taken as the obvious starting point for estimating surface adsorption potentials 
(ref. 15) . Application to new adatoms for which the interaction forces are unknown 
usually relies on atomic spectroscopic data to estimate atomic polarizabilities 
(ref. 16) . Although forms equivalent to the Kirkwood-Muller formula are used for 
dipole-quadrupole and quadrupole-quadrupole interaction terms (ref. 17) , the 
Kirkwood-Muller formula is still used to estimate the dispersion force constant. 

More recently, methods have been developed which allow calculations of inter- 
action potentials between closed-shell atoms (refs. 18 and 19). These calculations 
agree satisfactorily with known noble gas interactions and provide new results for 
noble gas interactions with alkali and halide ions. Since we are interested in thin- 
film growth on salts, which requires interactions with salt ions, these interactions 
seem to be a natural starting point for a systematic study of potential parameters 
to provide a data base for estimates of physical adsorption on salts. 

In the present report, we examine interatomic forces using semiempirical notions 
and derive relations for potential parameters in terms of basic atomic properties 
and repulsive core parameters. We require the formalism and the repulsive core 
parameters to describe diatomic interaction data, noble gas crystal data, and ionic 
crystal data. Such a globally consistent theory and data base should then provide a 
reliable basis for the estimation of potentials among many atom-atom and atom- 
surface systems. Because the potential parameters are related to fundamental atomic 
constants, the formation of a predictive capability is an inherent part of the 
present formalism. 


INTERATOMIC POTENTIAL 

The repulsive core was studied by Zener (ref. 9) , who shows that the quantum 
analog is 


4> R (r) ~ A exp(-br) 


( 1 ) 


where the exponential factor for two unlike atoms is found according to the combina- 
tion rule : 


AB 


- !< b * + b u> 


( 2 ) 


for atoms A and B. The short-range forces between noble gas atoms are also calcu- 
lated in the electron gas model by Waldman and Gordon (ref. 19) as well as by Gilbert 
and Wahl (ref. 20) using the Hartree-Pock methods. These investigations confirm the 
repulsive potential form given by Zener. As shown by Gordon and Kim (ref. 18) , 
repulsion is mainly made up of the total electronic kinetic energy as a function of 
nuclear separation and is closely related to the product of the atomic wave functions 
in the overlap region as assumed by Zener. Asymptotic wave functions have been found 
by Handler, Smith, and Silverstone (ref. 21) to be of the form 
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( 3 ) 



exp (-£ . r) 


where 


*i 
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and is the nuclear charge and is the number of electrons. Herein, is 

taken as a free parameter for each atomic species. The extension of Zener's repul- 
sive term is then 


<f> R (r) 


A r 


exp (-br) 


(5) 


with 

X = X A + X fi (6) 

and 

b = 5 a + C B (7) 

The are assumed to be closely related to the ionization potential and are to be 

determined. 

The long-range attractive- force potential is well approximated by induced 
electrostatic moments caused by the perturbing presence of the second atom. At 
distances large enough that the electron wave functions do not significantly overlap 
but not so large that the signal delay times between the atoms are unimportant, the 
long-range force potential may be expanded as 
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where the C L £ are related to the electric multipole polarizabilities (refs. 22 
and 23) given by 
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( 11 ) 
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where a L and a'^ are static polarizabilities, Z^ and Zg are atomic numbers, 

and N a and N B are the number of electrons in atoms A and B. These terms are 
equal to zero unless one of the atoms has a net charge imbalance (i.e., each expres- 
sion contains a monopole term) . The usual dispersion terms are given by (ref. 22) 
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where ( ioo) and a^(iu)) are the atomic dynamic polarizabilities. It has also 

been shown that a reasonable estimate of dispersion properties is given by the single 
frequency dispersion relation 
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for which equations (12a) to (12d) are approximated by 
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( 14c) 
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Hence, the attractive forces are known in terms of atomic spectral properties. We 
now require a formalism by which the required atomic dispersion parameters may be 
estimated. 


Dispersion Parameters for Hydrogen 

Dispersion relations for the hydrogen-like atom are exactly known (refs. 22 
and 23) . The static polarizability is 



2 ( 2L + 2) ! (L + 2) 
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Similarly, the asymptotic limit (to -*■ 00 ) of dynamic polarizability is given as 
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Using equations (15) and (16) we find the lowest energy poles satisfy 
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As is well known, the value of ( 0 ^ is taken from 
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(19) 


= 6.5 



Through use of the above relations , the dispersion force constants shown in table I 
are given to within 0.2 percent for the H-H interaction calculated elsewhere 
(ref. 22) . 


Dispersion Parameters for Two-Electron Systems 

The dipole spectral properties of atoms and dipole-dipole dispersion forces are 
well known for many atoms (refs. 23 to 25). However, we desire the means to estimate 
the dispersion forces on the basis of available data. The static dipole polariza- 
bility for helium is given by (refs. 26 and 27) 


a 2 He = 1.38 


( 20 ) 


whereas ^ = 1-46 (ref. 26) , yielding 


W, He =1.02 


( 21 ) 


The remaining a T and W T are estimated from (refs. 22 and 23) 


H A 
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( 22 ) 


and the frequency ratios of equations (17) and (18) . Values for Li and H are 
similarly obtained. The force coefficients in table I compare well with values found 
by Davison (ref. 23) . 


Dispersion Parameters for Noble Gases 

The dipole polarizability for noble gas atoms has been evaluated by Starkschall 
and Gordon (ref. 26) and values are shown in table II. The dipole dispersion energy 
pole is found with equation (14a) . Quadrupole and octapole polarizabilities are 
found in Doran (ref. 28) along with their dispersion force constants. It remains for 
us to estimate dispersion energy poles from the static multipole polarizabilities and 
their contributions to the dispersion forces . This is accomplished through use of 
equations (14a) to (14d) . Because some inconsistencies appear for the Xe data of 
Doran (ref. 28) , it is first necessary to observe that the spectral properties of Kr 
and Xe are similar before utilizing the value of the dipole energy pole of Xe. 


Dispersion Parameters for Ions 

The static dipole polarizabilities have been calculated from the numerical solu- 
tions of the Schrodinger equation (ref . 29) and are shown for halogen ions through 
chlorine and for alkali metal ions through potassium in table III along with the 



corresponding noble gas values . 
sequence have the same value of 


We assume that the atoms of the same isoeiectronic 
V, where 


V 



( 23 ) 


and I is the ionization potential. The parameter V is determined from the noble 
gas dispersion parameter and known values of I (ref. 30) . The dipole and quadru- 
pole polarizabilities of the alkali ions are taken from Dalgarno (ref. 27) . The 
remaining values for quadrupole polarizabilities for F - and Cl - are found in analogy 
to harmonic oscillator dispersion parameters to satisfy 


a 4 * 1.62 


a. 


W-, 


(24) 


and are assumed here for Br and I - as well. 

The dispersion forces have been evaluated in this paper with harmonic oscillator 
spectral distributions to obtain (ref. 31) 




= 210 < W ' 2/~1 

1,3 8 + 3U)^ 


(25) 


which is to be compared with the single pole dispersion (eq. (14c) ) 
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Approximate values of a g can be found by equating (25) and (26) , which results in 


1.56 (oo-^ + 0J 3 )a 2 
S 3“l 


(27) 


where the ionic values of t0 g are found from the corresponding noble gas values of 
the same isoeiectronic sequence as 
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( 28 ) 


The final values for the alkali and halide ions are shown in table III in comparison 
with the noble gas parameters. 



Ionic Electric Field 


As an atom comes close to an atomic ion, the Coulomb field experienced by the 
atom is substantially different from what it experiences from a point charge. To 
determine more accurately the electrostatic field we solve Poisson's equation for the 
ionic charge distribution as 


- — -or[r<t> (r ) ] = 4irp„ (r) 
r dr 2 


(29) 


where p g (r) is the electron distribution. The electron wave functions are taken 
as Slater's screen wave functions (ref. 32), which are approximately applicable to 
ions (ref. 33) . The screening factor c is given in table IV along with the corre- 
sponding principal quantum number n. The corresponding electrostatic potentials 
are 


c})^(r) = ~(cr + 1) exp(-2cr) + ^ 


(30) 


* 2 (r) 


(cr) J , (cr) 
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exp(-cr) + 


(31) 


and similarly for <j> n (r), where Q is the net charge, 
are approximated by scaling the Is potentials as 


The higher ordered potentials 


L (r) = -^/l + — \ exp (-2cr/n 2 ) + § 


(32) 


The corresponding electric field strength is 


E„(r) = 
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1 + 


2cr 


2cr 


exp(-2cr/n^) + 


(33) 


where Q = 1 for a positive ion and -1 for a negative ion. (Note that in these 
approximations the inner shell electrons are treated as negative point charges at 
the nucleus.) 


General Potential Form 

According to the calculations previously discussed, the general form for the 
interatomic potential is 
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where only A and b are as yet to be determined for each pair combination. The 
parameters of equation (34) are given by 


a 4 Qi 2 + «4 Qi 2 

i J i 

c 6 “ 2 + c l,l 

c 8 = C 1 , 2 + C 2 , 1 

C 10 = C 1 , 3 + C 3 , 1 + C 4 , 4 
The unknown potential parameters A and 
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b may be related through the following: 
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Equation (38) results in 
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(38) 
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whereas equation (39) results in 
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It is noted here that 


b 



(42) 


In the subsequent analysis, we search for values of A and b using the dispersion 
force and electric field parameters in tables I to IV. 

Force Parameter Combination Rules 

Just as the dispersion force parameters of atoms can be combined, so can the 
remaining force parameters referring to the core repulsion. The potential minimum 
occurs at the parameter value R for which the combination rule is given as 


R ij 



Rjj ) 


(43) 


Values for R^ are given in table V. 


<f> R (r) = A exp (-$) r^ exp 3 [1 


so that the combination rule for 3 is 


The core repulsion is written as 
(r/R) ] 


(44) 
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The A^j and 3-^ j are determined from an analysis of the Waldman and Gordon 
data (ref. 19) . We solve for (3^ by requiring that the potential evaluated at R^j 

gives the value -e ij and the corresponding derivative of the potential is zero, 
thus yielding 


3. . 

ri 
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(46) 
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The Waldman and Gordon data are then divided into a group involving neutral atoms 
only. In the case of i = j , we may solve for 6^ directly; for i ^ j, values 
are obtained by simultaneous iteration. The least-squares error value is selected 
from these results. The condition that i be a neutral atom and j a positive or 
a negative ion allows us to solve for 3jj directly. The values obtained are a 
function of and consequently we choose to minimize the sum of squared 

errors associated with all the Waldman and Gordon data. The data base is then 
expanded to include Herzberg's diatomic alkali-halide molecule data as compiled by 
Huheey (ref. 34) . Again, the values of are adjusted to minimize the sum of 

the squared errors of the enlarged data set. The next step in the parameter search 
is to include crystal data. 


NOBLE GAS SOLIDS 

The many-body noble gas system is described by the first-order Hamiltonian as 


H - + < 4 

i r j 

where cj)^j(r^j) is the two-body potential and is the kinetic energy operator. 

The cohesive energy E cq ^ is given as the sum of the zero point kinetic energy E z 
and the potential at the corresponding lattice site: 


E coh = % + S hj ( q> < 

J 

where rj is the distance to all the surrounding lattice sites. The zero point 
energy is related to the Debye temperature and is given by Glyde (ref. 35) . It is 
customary to take the sublimation energy L Q as the separation energy from a 
surface kink site. If we ignore the surface energy, L q is one-half the cohesive 
energy: 


L 


o 



(49) 


The experimental sublimation energy (ref. 36) and lattice spacings are used to 
further refine our values of and Rj_j_ for the noble gases. The two-body 

potential data are shown in table VI, and final values for noble gas potentials 
are shown in table VII. 

From the radii of the heavy noble gas atoms in table V, it is observed that 

R. 1 is proportional to \fl. From the expression of Zener (ref. 9) for the repulsive 

potential, one can expect that is nearly constant for at least the heavy noble 

gas atoms. The final values of R and £ for the noble gas systems are shown in 
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table VI with values for the Lennard-Jones potential, which is still used in many 
studies (e.g., ref. 16), and the range of results found in Waldman and Gordon 
(ref. 19). 


IONIC CRYSTALS 

The potential parameters are further refined by adding alkali-halide crystal 
data to the data base. Utilization of this added data is now described. 

The many-body crystal Hamiltonian (eq. (47)) is written as (ref. 37) 


H = £ K i + f £ 4>ij< r ij> 
i/ j 

where is the two-body potential appropriate for the i,j pair. The 

cohesive energy is given by 
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ID 


(50) 


where r — are taken as the equilibrium lattice position vectors and E z is the 
zero point energy of the crystal (ref. 38). The potentials appearing in equa- 
tion (50) are somewhat different from the results given by equation (34) . The 
reason for this difference is the induction force terms do not appear in equa- 
tion (50) since the electric fields at each lattice site are zero. The effective 
potential within the crystal at a given lattice site is thus 
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(51) 


where the induction terms of Cg are likewise zero. 

In the usual analysis of salt crystals, the potential is calculated at a given 
charge site and it is assumed that the sum is independent of the site charge 
(refs. 11, 37, and 39 to 44). This is accomplished by neglecting the repulsive terms 
for like ions (refs. 37 and 41 to 44) or by assuming the repulsive terms for like 

ions have functional forms that are identical, for example, A ++ = X = A + _ and 

b ++ = b = b + _ in equation (51) (refs. 11, 39, 40, 42, and 43). In reality, the 

crystals with large negative ions, such as Br~ and I - , or with unusually small 
positive ions, such as Li + , cannot be analyzed in this fashion because the negative 
ion repulsive cores are quite significant in determining the crystal size. For this 
reason, we use explicit crystal sums rather than one of the usual simplified formulas 
(refs. 11, 37, and 39 to 44). 
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The crystal lattice parameters are taken from the compilation of Fumi and Tosi 
(ref. 39) and the cohesive energy values are taken from the thermodynamic parameters 
(refs. 37, 39, 40, and 45). The atomic radii of the ions and the repulsive core 
parameters 3 -q are chosen as the best fit for the crystal data and Herzberg's 
diatomic molecular data as compiled by Huheey (ref. 34). In this process, an 
adequate representation of the noble gas and ion potentials of Waldman and Gordon 
(ref. 19) is maintained. (See table VIII.) The crystal data and data from the 
present model are shown in table IX. 


INTERACTIONS OF NOBLE GAS WITH CRYSTAL SURFACES 

Sufficient data now exist for the study of the interaction of noble gas atoms on 
noble gas solids and salt crystal surfaces. The present calculations are performed 
without taking into account the surface distortions related to the solid-surface 
energy. In this way the surface energy effects may be evaluated by performing the 
more exact calculations . 

The surface potentials for noble gas atoms on the noble gas crystal of the same 
element were evaluated for the three fee surfaces (100) , (110) , and (111) . In 

addition to evaluation of the potentials for the plane, two configurations were 
introduced on each of the planes to more accurately model a real surface. A half- 
plane sheet was added to study the effects of large aggregations (edge) , and a kink 
in the edge was also introduced and studied. To complete the information needed for 
kinetic studies, the saddle point energy between stable points on a flat surface was 
also evaluated. The results are shown in table X. For comparison, the two-body 
potential well depths £ are also given. 

It is shown in table X that the atoms of the noble gas crystal are most mobile 
on the close-packed (111) surface. The bonding at an edge is somewhat weak unless 
there is a kink site. Dimer formation at an edge is quite stable on the (111) 
surface. The (100) and (110) surfaces have more stable trapping sites with much 
less surface mobility than the (111) surface. Dimer formation on the (100) and (110) 
surfaces is quite stable. There is little difference between trapping at an edge 
and at a kink site on the (100) and (110) surfaces. 

Numerous calculations of the adsorption of noble gas atoms onto alkali-halide 
surfaces have been made over the years, each with its own variation in the inter- 
action potential (refs. 46 and 47). For example, table XI reveals the disparity in 
several of the more recently published theoretical calculations (refs. 48 and 49) 
for zero coverage adsorption potentials for the Ar/NaCl(100) system. Values are 
presented for four sites on the NaCl surface for the relaxed and unrelaxed state: 
a is the saddle point; b is the midpoint of the subcell edge; c is the cation; 
and d is the anion. Values shown in parentheses are potentials obtained from 
experimental heats of adsorption independent of surface heterogeneity and lateral 
interaction. These values are placed next to the calculated values of site c since 
it is the preferred site for Ar/NaCl (100) . The relationship between heat of 
adsorption AH q and the calculated potential p is given by 

-AH q = -c|> - E z + RT (52) 
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The variation in these calculated values is basically the result of any combination 
of (a) the selected potential forms and parameters, (b) neglect of the ions' finite 
size, (c) the number of Ar and alkali-halide ion pairs summed, and (d) the state of 
surface relaxation. Considering the values for an unrelaxed surface and the fact 
that summations for most of the papers examined were in excess of 2000 pairs, the 
discrepancies are still quite large. It is therefore likely that the major problem 
is associated with the selection of the specific form for the interatomic potentials 
and finite size corrections to the ion charge distribution (ref. 50) . Some improve- 
ment over past efforts may have been achieved for the work presented herein, since 
there is some overall consistency with published values (i.e., <7 percent error) in 
the global analysis of the paired interactions and of noble gas and alkali-halide 
crystal data. 

Unfortunately, there is little or no modern experimental data available with 
which these calculations can be compared. The work of Hayakawa (ref. 46) is perhaps 
the most thorough experimental work available, but it was conducted under very poor 
experimental conditions compared with present-day surface science techniques . For 
example, Hayakawa' s experiments appear to have been made at an ultimate pressure p 
of greater than 10“ ^ torr. At this pressure the surface of the alkali-halide 
crystals would have been completely covered by the chemical adsorption of residual 
gases such as CO, C0 2 , H 2 0, H 2 , and other adsorbates. Physical adsorption of the 
noble gases would have occurred on top of this intermediate layer and therefore 
would have substantially altered the true interaction potential. Furthermore, con- 
taminants from sample preparations such as polishing and exposure to atmosphere were 
not removed or even determined. In a more recent paper on Ar/NaCl (100) by Jackson 
and Davis (ref. 47) , some effort was made to improve the quality of the (100) 
crystals by double heat treatment, but in every other respect the approach was 
similar to Hayakawa' s. In table XII, the interaction potentials for an adatom of 
Ar on the ideal (100) surface of 16 alkali-halides are presented. 

Gas surface diffraction is another area of experimental research that provides 
potentials that can be compared with the calculations in this work. The bound state 
resonance spectra that are obtained from diffraction experiments can be related to 
some appropriate model potential with two or more free parameters. Adjusting these 
parameters to fit the resonance spectra then allows a determination of either the 
well depth or the interaction potential. This technique is very sensitive to the 
model used and can generate large disparities in well depths for the same experi- 
mental data. Calculations with a continuum model by Vidali, Cole, and Klein 
(ref. 51) for Ar/Ar(100) and Ar/Ar(lll) give values of 862.5 and 954.5 cal/mole, 
respectively, compared with our saddle point values of 1087.0 and 1090.0 cal/mole, 
respectively. Hoinkes (ref. 52) has shown that the partially summed continuum 
model yields well depths that are less than but near that of the saddle point for 
a Lennard-Jones (12-6) pair potential sum. Considering that the potential used in 
this work gives well depths that are 32 percent deeper than Lennard-Jones' values 
for Ar/Ar (see table VI), the agreement between our values and those of Vidali, 

Cole, and Klein would be improved had they used our interaction potential. Simi- 
larly, calculations are presented in table XIII for He on several ideal alkali- 
halide (100) surfaces at the three primary sites. Available well depths from 
scattering experiments are given in parentheses for comparison (refs. 51 to 53) . 

As previously stated, these values are sensitive to the model used for extracting 
the well depth and would probably scale in the direction which substantially 
improves agreement if our model were used. 
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CONCLUDING REMARKS 


The development of a theory and associated data base for the physical inter- 
action of two atoms has been achieved. The base draws heavily on many years of study 
of atomic dispersions, two-body potentials, noble gas crystals, and salt crystals. 

The data base shows an overall consistency with crystal data to within several per- 
cent. Hence, a reliable data base for surface physics studies is herein provided. 


NASA Langley Research Center 
Hampton, VA 23665-5225 
March 3, 1986 
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SYMBOLS 


atom of diatomic pair 
lattice parameter, atomic units 
Bohr's radius, 0.592 x 10~® cm 
atom of diatomic pair 
screening parameter 

zero point energy of lattice, cal/mol 
atomic ionization potential, atomic units 
sublimation energy, cal/mol 
valence shell principal quantum number 

equilibrium separation of two-body potential, atomic units 

potential energy of crystal, cal/mol 

distance above surface, atomic units 

multipole polarizability (2 L -pole) , (atomic units) 

repulsive core parameter 

well depth, atomic units 

ratio of co^/I 

diatomic potential energy at separation distance r, atomic units 
pole location of the dynamic 2 L -polarizability , atomic units 



TABLE I.- DISPERSION PARAMETERS AND FORCE CONSTANTS 
[values in atomic units] 


Parameter 


Atomic dispersion parameters 


4.5 

0.428 

15.0 

0.574 

131.3 

0.651 


0.281 

174 

1.38 

1.71 

0.077 

1.02 

0.234 

5 130 

2.44 

2.30 

0.103 

1.37 

0.512 

33 800 

11.2 

2.61 

0.117 

1.50 


Hydrogen atom forces 


0.192 

3.10 

0.114 

4.16 

0.179 

4.72 


C l,l 

6.50 

0.65 

76.7 

2.81 

0 . 

487 

C l,2 

62.1 

1.43 

7 187 

13.4 

0 . 

747 

C 2,l 

62.1 

6.8 

664 

28.5 

5. 

23 

C 1 / 3 

640.8 

3.58 

58 696 

70.5 

1 . 

33 

C 3 , 1 

640.8 

7.3 

6 607 

30.2 

57 

.0 

C 3 , 2 

1130 

28.3 

117 599 

259 

15 

.1 



Helium atom 

forces 




0.372 

25.8 

1.46 

0.86 

2 483 

7.38 

1.96 

116 

7.38 

2.19 

20 562 

39.4 

10.6 

599 

39.4 

8.6 

20 984 

71.4 


0.305 

0.483 

1.67 

0.87 

9.13 

5.02 





















TABLE III.- DISPERSION PARAMETERS FOR CLOSED-SHELL IONS AND NOBLE GASES 


[values in atomic units] 


Atom 

a 2 

I 

V 

“1 

a 4 

w 2 

a 8 

W 3 

F - 

6.75 

0.11 

1.53 

0.168 

57.5 

0.182 

526 

0.41 

Ne 

2.67 

0.793 

1.53 

1.21 

6.42 

1.09 

30.4 

2.96 

Na + 

1.55 

1.68 

1.53 

2.57 

1.57 

2.32 

3.64 

6.27 

Cl" 

22.0 

0.10 

1.22 

0.122 

325.9 

0.126 

3585 

0.22 

Ar 

11.1 

0.579 

1.22 

0.701 

- 48.2 

0.504 

449 . 7 

1.30 

K + 

6.38 

1.13 

1.22 

1.38 

13.6 

0.97 

56.9 

2.54 

Br ~ 

27.4 

0.09 

1.18 

0.106 

418 

0.116 

5242 

0.29 

Kr 

16.8 

0.515 

1.18 

0.605 

78.8 

0.455 

681.4 

1.67 

Rb + 

6.75 

1.02 

1.18 

1.204 

~15 

0.91 

69.8 

3.31 

1 “ 

41.9 

0.12 

1.20 

0.144 

471 

0.111 

4287 

0.40 

Xe 

27.1 

0.446 

1.20 

0.536 

128.3 

0.403 

1105 

1.48 

Cs + 

18.2 

0.86 

1.20 

1.03 

~49 

0.77 

257 

2.85 











Atom 

or 

ion 

R, atomic 
units 

He 

5.8200 

Ne 

6.0109 

Ar 

7.0632 

Kr 

7.5000 

Xe 

8.1515 

Li + 

2.2130 

Na + 

3.4100 

K + 

4.5940 

Rb + 

5.0585 

F“ 

4.5280 

Cl" 

6.9276 

Br - 

7.1397 

1“ 

7.7199 


3, atomic 
units 

14.100 

15.523 

15.104 

15.883 

15.268 

5.097 

8.575 

10.585 

13.050 

13.251 

15.758 

19.692 

30.776 







Lennard- Jones 
values 

Present values 

Reference 19 values 

System 

R, 

atomic 

units 

e, 

atomic 

units 

R, 

atomic 

units 

atomic 

units 

R, 

atomic 

units 

e, 

atomic 

units 


He-He 

5.834 

2.03 

X 

10" 5 

5.820 

2.80 

X 

10“ 5 

5.59-5.82 

3.48-3.58 

X 

10“ 5 

He-Ne 

5.921 

4.15 

X 

10“ 5 

5.915 

5.62 

X 

10“ 5 

6.01-6.07 

4.52-5.48 

X 

10“ 5 

He-Ar 

6.517 

7.53 

X 

10“ 5 

6.442 

1.07 

X 

10“ 4 

6.54-6.90 

6.45-11.6 

X 

10“ 5 

He-Kr 

6.740 

88.7 

X 

10“ 5 

6.660 

1.24 

X 

10“ 4 

6.67-7.09 

7.83-16.4 

X 

10“ 5 

He-Xe 

7.067 

9.41 

X 

10“ 5 

6.986 

1.34 

X 

10“ 4 

7.09-7.84 

7.98-17.9 

X 

10“ 5 

Ne-Ne 

6.008 

8.59 

X 

10“ 5 

6.011 

1.13 

X 

10~ 4 

5.88-6.05 

1.12-1.35 

X 

10“ 4 

Ne-Ar 

6.604 

1.51 

X 

10“ 4 

6.537 

2.13 

X 

10“ 4 

6.48-6.77 

1.84-2.27 

X 

10“ 4 

Ne-Kr 

6.827 

1.78 

X 

10~ 4 

6.755 

2.45 

X 

10“ 4 

6.77-7.03 

2.20-2.36 

X 

10“ 4 

Ne-Xe 

7.152 

1.88 

X 

10“ 4 

7.081 

2.65 

X 

10" 4 

7.09-7.56 

2.11-2.38 

X 

10“ 4 

Ar-Ar 

7.200 

3.00 

X 

10“ 4 

7.063 

4.39 

X 

10“ 4 

7.11-7.11 

4.50-4.61 

X 

10“ 4 

Ar-Kr 

7.423 

3.61 

X 

10“ 4 

7.282 

5.18 

X 

10“ 4 

7.29-7.33 

5.18-5.78 

X 

10“ 4 

Ar-Xe 

7.748 

4.00 

X 

10“ 4 

7.607 

5.77 

X 

10“ 4 

7.65-7.75 

6 . 01-6 . 68 

X 

10“ 4 

Kr-Kr 

7.646 

4.38 

X 

10“ 4 

7.500 

6.17 

X 

10“ 4 

7.50-7.58 

6.38-7.00 

X 

10 4 

Kr-Xe 

7.971 

4.91 

X 

10“ 4 

7.826 

6.95 

X 

10“ 4 

7.81-7.90 

7.25-8.37 

X 

10“ 4 

Xe-Xe 

8.295 

5.63 

X 

10 -4 

8.152 

7.93 

X 

10“ 4 

7.99-8.26 

8.24-11.0 

X 

10“ 4 























TABLE VII.- NOBLE GAS PAIR INTERACTION PARAMETERS AND SOLID PARAMETERS 
FROM PRESENT MODEL AND EXPERIMENTS 


Gets 

type 

cal/mol 

(a) 

R, 

atomic 

units 

(a) 

a o' 

atomic 

units 

(a) 

a o ' 

atomic 

units 

(b) 

E Z ' 

cal/mol 

(a) 

L o' 

cal/mol 

(a) 

E o ' 

cal/mol 

(b) 

Ne 

70.7 

6.011 

8.438 

8.432 

149 

468 

473+34 

Ar 

274.6 

7.063 

9.866 

10.038 

208 

2010 

1864+34 

Kr 

386.0 

7.500 

10.502 

10.676 

160 

2864 

2659+12 

Xe 

496.1 

8.152 

11.367 

11.586 

123 

3798 

3749+58 


a Present model. 
^Experiments . 
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TABLE VIII.- DIATOMIC INTERACTION POTENTIAL PARAMETERS £ AB AND r AB 
FROM PRESENT MODEL AND LITERATURE 

[All quantities are in atomic units] 


A 

B 

e AB from 
present 
model 

R AB f rom 
present 

model 

£ AB f rom 
literature 

R ab from 
literature 

Source 

Li + 

He 

1.58 

X 

10" 3 

4.017 

1 

.72-2. 

76 

x 10 -3 

3.666-4.195 

Reference 

19 


Ne 

3.02 

X 

10" 3 

4.112 


4.52 

X 

10" 3 

3.704 

Reference 

19 


Ar 

8.00 

X 

10" 3 

4.638 


1.13 

X 

10" 2 

4.176 

Reference 

19 


Kr 

10.3 

X 

10" 3 

4.857 


1.43 

X 

10" 2 

4.365 

Reference 

19 


Xe 

12.7 

X 

10" 3 

5.182 


1.67 

X 

10" 2 

4.743 

Reference 

19 


F 

2.85 

X 

10 , 

3.371 


2.87 

X 

10" 1 

2.923 

Reference 

34 


Cl" 

2.22 

X 

10 t 

4.570 


2.40 

X 

10" 1 

3.817 

Reference 

34 


Br 

2.19 

X 

10 1 

4.676 


2.30 

X 

10' 1 

4.101 

Reference 

34 


I" 

2.06 

X 

10" 1 

4.966 


2.13 

X 

10" 1 

4.520 

Reference 

34 

Na + 

He 

1.16 

X 

10" 3 

4.615 

1 

.29-1. 

47 

x 10" 3 

4.441-4.573 

Reference 

19 


Ne 

2.21 

X 

10" 3 

4.710 


2.42 

X 

10" 3 

4.592 

Reference 

19 


Ar 

5.86 

X 

10" 3 

5.237 


5.92 

X 

10" 3 

5.140 

Reference 

19 


Kr 

7.57 

X 

10" 3 

5.455 


7.43 

X 

10" 3 

5.367 

Reference 

19 


Xe 

9.45 

X 

10 3 

5.781 


8.86 

X 

10" 3 

5.764 

Reference 

19 


F" 

2.44 

X 

10" 1 

3.969 


2.41 

X 

10" 1 

3.477 

Reference 

34 


Cl" 

1.94 

X 

10" 1 

5.169 


2.08 

X 

10" 1 

4.462 

Reference 

34 


Br" 

1.93 

X 

10" 1 

5.275 


1.99 

X 

10" 1 

4.728 

Reference 

34 


1“ 

1.80 

X 

10" 1 

5.565 


1.87 

X 

10" 1 

5.125 

Reference 

34 

K + 

He 

8.98 

X 

10" 4 

5.207 


8.13 

X 

10" 4 

5.442 

Reference 

19 


Ne 

1.72 

X 

10" 3 

5.302 


1.43 

X 

10" 3 

5.537 

Reference 

19 


Ar 

4.43 

X 

10" 3 

5.829 

4 

.08-5. 

14 

X 10" 3 

5.442-5.896 

Reference 

19 


Kr 

5.69 

X 

10" 3 

6.047 


5.40 

X 

10" 3 

6.066 

Reference 

19 


Xe 

7.11 

X 

10" 3 

6.373 


6.88 

X 

10" 3 

6.350 

Reference 

19 


F" 

2.15 

X 

10" 1 

4.561 


2.16 

X 

10" 1 

4.025 

Reference 

34 


Cl" 

1.73 

X 

10" 1 

5.761 


1.84 

X 

10" 1 

5.044 

Reference 

34 


Br - 

1.71 

X 

10" 1 

5.867 


1.76 

X 

10" 1 

5.331 

Reference 

34 


I" 

1.61 

X 

— 1 
1 

o 
1 — 1 

6.157 


1.66 

X 

10" 1 

5.760 

Reference 

34 
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A 

B 

from 

present 

model 

R A B from 
present 
model 

e AB from 
literature 

Rab from 
literature 

Source 

Rb + 

He 

8.58 x io~ 4 

5.439 

7.54 x io -4 

5.726 

Reference 19 


Ne 

1.64 x io -3 

5.535 

1.25 x io -3 

5.915 

Reference 19 


Ar 

4.21 x io -3 

6.061 

3.63 x io -3 

6.236 

Reference 19 


Kr 

5.39 x io -3 

6.279 

4.38-4.67 x io -3 

6.312-6.406 

Reference 19 


Xe 

6.74 x io -3 

6.605 

6.03 x io -3 

6.609 

Reference 19 


F~ 

2.08 x io _1 

4.793 

2.10 x IO" 1 

4.282 

Reference 34 


Cl" 

1.68 x 10 -1 

5.993 

1.86 x icT 1 

5.267 

Reference 34 


Br~ 

1.66 x io -1 

6.099 

1.72 x io -1 

5.565 

Reference 34 


I" 

1.56 x 10 -1 

6.389 

1.62 x 10 -1 

6.004 

Reference 34 

P" 

Xe 

7.63 x io -3 

6.340 

8.38 x 10 -3 

5.820 

Reference 19 

Cl" 

Xe 

4.02 x io“ 2 

7.540 

6.36 x io -3 

6.595 

Reference 19 

1“ 

He 

2.73 x io -4 

6.770 

3.20 x io~ 4 

7.578 

Reference 19 


Ne 

5.17 x io -4 

6.865 

3.09 x io -4 

8.485 

Reference 19 


Ar 

1.55 x io -3 

7.392 

2.02 x io -3 

7.597 

Reference 19 


Kr 

2.05 x lo“ 3 

7.610 

2.94 x io -3 

7.597 

Reference 19 


Xe 

2.74 x IO -3 

7.936 

5.29 x io“ 3 

7.276 

Reference 19 

Br - 

He 

3.69 x io -4 

6.480 

3.57 x io -4 

7.105 

Reference 19 


Ne 

6.96 x lo“ 4 

6.575 

4.08 x io" 4 

7.861 

Reference 19 


Ar 

2.08 x IO -3 

7.101 

2.32 x io“ 3 

7.124 

Reference 19 


Kr 

2.76 x io -3 

7.320 

3.42 x io -3 

7.105 

Reference 19 


Xe 

3.67 x lcT 3 

7.646 

5.84 x io -3 

6.879 

Reference 19 














TABLE IX.- IONIC CRYSTAL PARAMETERS AND EXPERIMENTAL VALUES 


Crystal 

a, 

atomic 

units 

u, 

kcal/mol 

E z , 

kcal/mol 

Experimental 

a, 

atomic ' 
units 

Experimental 

u, 

kcal/mol 

Source of 
experimental data 

Li + F“ 

3.932 

258 

2.81 

3.802 

246±3 

References 

39,40,45 

Li + Cl“ 

5.251 

206 

1.43 

4.858 

202±3 

References 

37,39,40,45 

Li + Br“ 

5.280 

205 

1.51 

5.193 

192±4 

References 

37,39,40,45 

Li + I“ 

5.518 

187 

.85 

5.669 

178±6 

References 

37,39,40,45 

Na + F~ 

4.491 

225 

1.92 

4.372 

218±3 • 

References 

39,40,45 

Na + Cl“ 

5.729 

184 

1.28 

5.323 

18613 

References 

37,39,40,45 

Na + Br~ 

5.767 

181 

.95 

5.641 

17714 

References 

37,39,40,45 

Na + I~ 

5.962 

165 

.79 

6.111 

166+3 

References 

37,39,40,45 

K + F~ 

5.145 

197 

1.42 

5.093 

193+3 

References 

39,40,45 

K + C1“ 

6.374 

162 

.81 

5.939 

168+3 

References 

37,39,40,45 

K + Br~ 

6.398 

161 

1.03 

6.228 

162+2 

References 

37,39,40,45 

K + I~ 

6.623 

148 

.67 

6.670 

15212 

References 

37,39,40,45 

Rb + F~ 

5.343 

194 

1.13 

5.320 

184+2 

References 

39,40,45 

Rb + Cl~ 

6.591 

159 

.82 

6.198 

163+2 

References 

37,39,40,45 

Rb + Br~ 

6.627 

157 

.63 

6.493 

15613 

References 

37,39,40,45 

Rb + I~ 

6.856 

143 

.52 

6.930 

14813 

References 

37,39,40,45 
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TABLE X.- SURFACE POTENTIALS FOR NOBLE GAS ATOMS ON NOBLE GAS CRYSTALS 


[values in parentheses are multiples of e] 


Face 

Surface 

Potential for 
Ne (e = 70.7 ca 1/mol ) , 
cal/mol 

Potential for 
Ar (£ = 274.6 cal/mol) , 
cal/mol 

Potential for 
Kr (£ = 386.0 cal/mol) , 
cal/mol 

Potential for 
Xe (£ = 496.1 cal/mol), 
cal/mol 

(111) 

Flat 

Edge 

Kink 

Saddle 

316 (4.47) 
475 (6.72) 
547 (7.74) 
285 (4.03) 

1198 (4.47) 
1807 (6.38) 
2099 (7.60) 
1090 (3.97) 

1664 (4.31) 
2517 (6.52) 
2910 (7.54) 
1504 (3.90) 

2193 (4.42) 
3299 (6.65) 
3806 (7.67) 
1998 (4.03) 

(100) 

Flat 

Edge 

Kink 

Saddle 

372 (5.26) 
533 (7.54) 
547 (7.74) 
280 (3.96) 

1416 (5.15) 
2039 (7.43) 
2088 (7.60) 
1087 (3.96) 

1973 (5.11) 
2844 (7.37) 
2911 (7.54) 
1495 (3.87) 

2584 (5.21) 
3716 (7.49) 
3807 (7.67) 
1987 (4.01) 

(110) 

Flat 

Edge 

Kink 

Saddle 

452 (6.39) 
533 (7.54) 
547 (7.74) 
288 (4.07) 

1733 (6.31) 
2046 (7.45) 
2099 (7.60) 
1103 (4.02) 

2396 (6.21) 
2853 (7.39) 
2910 (7.54) 
1522 (3.94) 

3163 (6.38) 
3731 (7.52) 
3806 (7.67) 
2029 (4.09) 




TABLE XI.- SURFACE POTENTIALS FOR Ar/NaCl (100) FOR RELAXED 
AND UNRELAXED SURFACES 


[Values in parentheses are potentials obtained from experimental 
isosteric heats taken from these works] 


Site 

(a) 

Relaxed 

state 

Unrelaxed state 

Source 

cal/mol 

z/a 

< t > » 

cal/mol 

z/a 

a 

2006 

1.16 



Reference 48 

b 

1790 

1.22 




c 

2420 

1.04 




d 

1481 

1.32 




a 

1066 

1.13 

1023 

1.12 

Reference 49 

b 






c 

1109 

1.11 

801 

1.71 


d 

751 

1.32 

780 

1.28 


a 



1768 

1.17 

Reference 46 

b 



1514 

1.25 


c 



1876 (1619) 

1.17 


d 



1401 

1.32 


a 



1119 

1.19 

This work 

b 



946 

1.26 


c 



1484 

1.10 


d 



734 

1.38 


a 



1098 

1.18 

Reference 47 

b 



962 

1.24 


c 



1096 (1749) 

1.17 


d 



819 

1.34 



a Sites: a - saddle point; b - midpoint of subcell edge; c - cation; 

d - anion. 
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TABLE XII. - ZERO COVERAGE ADSORPTION POTENTIALS FOR Ar 
ON (100) SURFACE OF 16 ALKALI -HALIDES 


Alkali-halide 

Cation 

<tb 

ca 1/mol 

Cation 

z/a 

Anion 

<t>, 

cal/mol 

Anion 

z/a 

Saddle 

point 

4>, 

cal/mol 

Saddle 

point 

z/a 

LiF 

2270 

1.438 

1757 

1.6 

2132 

1.488 

LiCl 

1536 

1.2 

644 

1.413 

1107 

1.35 

LiBr 

2227 

1 . 025 

821.8 

1.4 

1310 

1.2 

Li I 

5790 

.838 

1227 

1.288 

2417 

1.063 

NaF 

1735 

1.313 

1518 

1.4 

1971 

1.263 

NaCl 

148.4 

1.1 

732.8 

1.375 

1119 

1.188 

NaBr 

1856 

.988 

763.7 

1.288 

1335 

1.075 

Nal 

3766 

.85 

1054 

1.213 

2268 

.95 

KF 

1376 

1.225 

1589 

1.188 

2033 

1.063 

KC1 

1314 

1.05 

774 

1.225 

1342 

1.013 

KBr 

1429 

.988 

799.2 

1.163 

1569 

.925 

RbF 

1144 

1.225 

1438 

1.15 

1775 

1.038 

RbCl 

1128 

1.039 

708.8 

1.175 

1279 

.95 

RbBr 

1219 

.975 

737.9 

1.125 

1523 

.863 

Rbl 

1958 

.888 

903.4 

1.063 

2300 

.775 


TABLE XIII.- INTERACTION POTENTIALS FOR He ADSORBED ON SEVERAL 
UNRELAXED ALKALI-HALIDE SURFACES 

[Well depths from scattering experiments given in 
parentheses for comparison] 


Primary 

site 

(j), cal/mol, for He adsorbed on - 

Source 

LF 

NaF 

NaCl 

Cation 

255.3 

236.4 

260.9 


Anion 

199.4 

210.8 

111.7 


Saddle point 

253.5 

293.4 

184.6 



(185.2) 

(162.2) 


Reference 51 


*(197.8) 

*(166.7) 


Reference 52 




(163.3) 

Reference 53 


^Represents average of numerous experiments. 
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